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Abstract 

Monte  Carlo  simulations  are  reported  for  the  solvent  primitive  model  (SPM)  of  the 
electrical  double  layer  at  electrically  charged  and  neutral  surfaces.  Both  solvent  and  ions 
aure  modeled  as  hard  spheres  with  interactions  governed  by  hard  sphere  and  Coulomb 
interactions.  Density  profiles  of  solvent  and  tons  and  electrostatic  potential  profiles  are 
presented  for  1:1  electrolyte  concentrations  of  1  and  2M  at  300K.  Comparison  of  results  at 
charged  and  neutral  walls  indicates  that  the  density  distribution  of  a  dense  solvent  near  a 
wall  induces  significant  layering  of  ions  and  decrease  of  potentials  in  the  double  layer.  This 
layering  effect  cannot  be  captiued  in  simulations  of  the  conventional  primitive  model  (PM) 
of  the  double  layer.  The  solvent  induced  steric  ordering  of  electrolyte  ions  at  a  neutral 
wall  can  create  a  space  charge  layer. 
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I.  Introduction 


Separation  of  charge  in  the  vicinity  of  a  surface  that  is  either  charged  or  neutral  occurs 
in  a  variety  of  phenomena  in  electrochemical,  biological,  tribolc^cal,  and  colloidal  sciences. 
For  low  surface  charge  and  dilute  electrolyte  solutions,  the  classical  theory  of  Gouy  (ij  and 
Chapman  [2],  and  later  modified  by  Stem  (GCS)  [3]  gives  predictions  of  the  interfacial  ion 
and  potential  distributions  that  are  in  excellent  agreement  with  experiment  and  computer 
simulations.  More  recent  theoretical  and  computational  studies  have  demonstrated  that 
the  detailed  stmeture  of  the  solvent  and  ion  are  fundamental  to  accurately  describe  the 
interfacial  double  layer  properties  for  highly  coupled  systems  at  high  bulk  ion  densities. 
Among  the  models  of  the  electrical  double  layer,  the  most  widely  used  is  the  primitive 
model  (PM),  in  which  ions  are  represented  as  point  charges  that  are  embedded  in  hard 
spheres,  the  solvent  is  modeled  as  an  isotropic  dielectric  continuum,  and  the  surface  mod¬ 
eled  as  a  hard  wall  with  uniform  surface  charge  density.  Vigorous  investigation  of  the  PM 
by  computer  simulations  [4*11],  integral  equation  methods,  and  functional  theories  [12*19] 
diiring  the  last  decade  have  revealed  phenomena  of  the  electrical  double  layer  that  are  not 
captured  by  either  the  GC  or  other  theories.  For  high  surface  charge  densities  or  high 
electrolyte  concentrations,  the  ion  denrity  profiles  in  the  PM  predict  a  highly  organized 
layering  of  ions  at  the  charged  surface  than  either  the  GC  or  GCS.  In  particxilar,  for  1:1 
electrolytes,  the  conterion  profile  shows  three  distinct  layers.  For  both  1:1  a^nd  2:2  elec* 
trolytes  at  high  concentrations,  there  is  charge  inversion  with  coion  densities  exceeding  the 
counterion  density  in  the  second  layer.  Neither  of  these  phenomena  are  predicted  by  the 
GCS  model,  although  the  probable  occiirrence  of  both  were  discussed  by  Grahame  [20]  in 
the  interpretation  of  the  interfacial  tension  measurements  of  Hg.  The  obvious  disadvan¬ 
tage  of  the  PM  is  that  the  interfacial  structure  of  the  solvent  is  minimized,  although  it 
is  well  known  that  ion  solvation  and  dipole  orientation  at  the  surface,  both  cS.  which  are 
critical  in  determining  the  potential  and  ion  distributions,  are  correlated  with  the  molecu¬ 
lar  structure  of  the  solvent.  The  development  of  realistic  models  of  water  and  the  civilized 
primitive  model  [21],  in  which  the  solvent  is  treated  as  a  dipole  embedded  in  a  hard  sphere, 
allows  the  effects  of  the  solvent  structure  on  the  double  layer  structure  to  be  studied.  The 
complexity  of  realistic  models  of  water  limits  studies  of  these  systems  to  computational 
simulations;  however,  very  recent  investigations  of  an  electrified  PilH-xO  interface  [22]  (no 
ions  present)  demonstrates  the  importance  of  molecular  solvent  structure  on  the  interface 
potential  distribution. 
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A  recent  theoretical  study  of  the  civilized  primitive  model  of  the  electrical  double  layer 
by  Russier  et  al.  [23]  notes  the  importance  of  steric  effects  imposed  by  the  solvent  on  the 
ion  densities  expected  at  a  hard  wall.  These  authors  showed  that  the  closest  approach  of  an 
ion  to  a  neutral  wall  is  determined  by  electrostatic  ion  dipole  interactions  between  ion  and 
solvent  which  is  proportional  to  the  inverse  square  of  the  distance  between  the  interacting 
hard  spheres  and  the  solvent  diameter.  The  results  of  Russier  et  al.  are  quantitatively 
different  from  any  previous  results  fotmd  in  the  PM,  in  that  the  closest  approach  of  the 
ion  does  not  neccesarily  correspond  to  intimate  contact  between  the  wall  and  ion. 

In  this  work,  we  present  Monte  Carlo  simulations  of  a  solvent  primitive  model  (SPM) 
of  the  electrical  double  layer  in  which  the  effects  due  to  the  finite  size  of  solvent  molecules 
are  included.  Ion  density  and  potential  profiles  are  compared  for  electrolytes  in  contact 
with  charged  and  neutral  wails  for  several  sirrface  charge  densities  and  electrolyte  concen¬ 
trations.  The  key  findings  of  our  present  work  is  that  the  finite  size  of  the  solvent  results 
in  highly  ordered  layering  of  ions  at  both  charged  and  neutral  walls.  For  electrolytes  com¬ 
prised  of  different  size  anions  and  cations,  the  solvent  induced  structure  results  in  a  space 
charged  layer  about  two  molecules  thick  at  a  neutral  wall. 


II.  Model  System 


The  systems  simulated  in  tins  work  are  mixtures  of  1:1  primitive  electrolytes  and 
hard  sphere  solvent  molecules  contained  between  two  infinite  planar  sxirfaces.  One  of  the 
surfaces  (at  z  =  0)  is  uniformly  charged  and  the  other  (at  z  =  H)w  neutral.  The  charged 
side  of  the  system  is  shown  schematically  in  Figure  1  for  solvent  molectiles  and  ions  of 
diameter  d.  In  the  computations,  the  systems  are  represented  as  a  periodic  array  of  closed 
rectangular  boxes  of  height  H  and  a  square  base  of  length  L  on  a  side.  The  surface  charges 
on  the  wall  at  z  s  0  is  balanced  by  the  excess  coimterions  in  the  solution.  and  iV_ 
and  q+  and  are  the  ntimber  and  charge  of  cations  and  anions  in  the  enclosed  system. 
Electroneutrality  requires 

|Ar_9_ -iV+^+l  =  (rA  ,  (1) 


where  cr  is  the  surface  charge  density  and  A  =  L^.  The  interaction  potential  between 


particles  i  and  j  is 


oo  , 
.im  , 

tln-rj  I  ’ 


nj  <  (di  +  dj)f2 
rij  >  (di  +  dj)/2  , 


(2a) 
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and  between  particle  i  and  the  wail  is, 


oo, 

--2ireqiZilt , 


Zi  <  di/2 
Zi  >  di/2  , 


(26) 


where  zi  is  the  distance  between  the  wall  and  particle  i.  In  equations  (2a,b),  di  is  the 
diameter  of  particle  i ,  g,  equals  the  ion  valence  times  the  unit  electronic  charge  t.  t  is  the 
dielectric  constant  (taken  to  be  78.5  ,  the  value  of  btilk  liquid  water,  in  the  simulation). 
Of  coxirse,  if  either  i  or  j  is  a  neutral  particle  the  Cotilomb  potential  in  (2a)  or  (2b)  is  zero. 

In  the  simulations  the  system  is  periodically  extended  in  the  x  and  y  directions.  The 
resulting  long  range  Coulomb  potential  between  ions  and  their  periodic  replicas  is  given 
by  the  sum  obtained  by  [7,24] 


OO 


=  4  cos(/0  +  2irm)*  +  C*)  -  ln(co8hC  -  coe rf). 


(3) 


with  the  notation  that  ^  =  2ir(rj  —  Xj)/L,  n  =  2fr(yi  —  yj)/L,  and  C  =  2x(zi  —  Zj)/L. 
Ko{x)  is  the  modified  Bessel  function  which  (^creases  to  zero  very  quickly  as  x  increases. 


III.  Simulation 


In  order  to  run  the  simtilation  corresponding  to  fixed  bulk  salt  concentrations,  we  first 
run  the  grand  canonical  Monte  Carlo  (GCMC)  simulations  for  the  primitive  electrolyte 
systems  in  the  following  way  [4-5].  The  addition  and  deletion  are  attempted  with  equal 
probability  of  about  10  percent  and  the  moves  are  attempted  at  80  percent  probability. 
At  each  addition  or  deletion  a  neutral  pair  of  cations  and  anions  are  added  or  deleted.  If 
the  probability  of  accepting  a  trial  step  firom  state  i  to  state  j  is  /t^,  then 


Nt\Nr\ 

Nf\Nj\ 


where  Nf  and  Nf  are  the  number  of  cations  and  anions  in  states  i  and  j  and 


(4) 


B  =!  2ln7±  +ln(n+T4)  +  ln(n-.VL)  , 


(5) 


with  7:^  the  mean  ionic  activity  and  nj  is  the  number  density  of  particle  i  with  =  Ni/Vy. 
Vi  is  the  volume  accessible  to  particle  i.  The  Markov  chain  is  set  up  according  to 


fij  =  m*n(l,  fij/fji)  additiariy 
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(6) 


fji  =  mtn(l,  fjiffij)  deletion, 

fij  sz  min[l,  exp(—0Uij)]  move, 

where  $  =  \/kT  with  k  is  the  Boltzmann  constant  and  T  is  the  temperature.  In  the 
simulations,  T  is  fixed  at  300  Kelvin.  The  ionic  activity  coefficients  ln7i;  =s  —0.127  and 
0.271  [5]  are  used  in  the  GCMC  simulations  respectively  for  systems  (o  —  c)  listed  in  Table 
I.  For  systems  d  and  e,  there  are  no  previously  available  data  for  ln7^.  Using  equation  (4), 
for  a  1  molar  1:1  electrolyte,  we  obtained  the  ionic  activities  from  bulk  electrolyte  GCMC 
simulations  to  be  ln7i;  —0.318  and  —0.44  for  systems  d  and  e,  r^pectively. 

For  the  confined  systems,  the  average  number  of  cations  and  anions  corresponding 
to  the  right  bulk  concentrations  of  the  systems  are  obtained  from  GCMC  simulations. 
Subsequently,  simulations  of  the  SPM  model  of  the  double  layer  were  carried  out  in  the 
canonical  Monte  Carlo  (CMC)  [7]  simulations  using  systems  (a  —  e).  The  total  number  of 
solvent  molecules  and  ions  was  fixed  at 

(N^  +  N^+  NQ)df^/HL‘*  -  0.7,  (7) 

where  and  iV_  are  the  average  number  of  free  ions  in  the  confined  PM  fluid.  This  total 
number  density  corresponds  roughly  to  an  aqueous  solution  at  300K.  At  such  high  density 
a  grand  canonical  Monte  Carlo  simulation  is  very  hard  to  carry  out  due  to  the  difficulties 
in  inserting  particles.  The  simulation  process  can  be  divided  into  four  steps  : 

(1)  GCMC  for  bulk  =>  ln7±  and  bulk  concentration  C 

(2)  GCMC  for  confined  system  with  ln7±  and  C  =►  Average  -f 

(3)  Add  solvent  =►  SPM,  nd*  ~  0.7 

(4)  CMC  for  SPM 

When  the  electrolyte  concentration  is  1  molar,  the  box  has  a  height  oi  H  ss  I2d  and 
base  of  £  =  4.4d;  reduced  surface  charge  densities  of  9*  =  o<P/e  =s  0.42  and  0.7,  with 
d  =  4.25  A,  were  itsed  in  the  simulations.  For  an  electrolyte  concentration  of  2  molar, 
H  =s  14.14d  and  L  s  4.5d,  and  o*  =  0.40.  These  parameters  as  well  as  the  solvent  and 
ion  diameters  are  listed  in  Table  I.  The  long  range  Coulomb  interaction  is  handled  by  the 
summation  technique  given  in  section  11.  The  probability  of  acceptanw  used  in  the  SPM 
simulations  is  10  percent.  During  che  run,  the  systems  were  equilibrated  for  10^  steps  and 
then  was  run  (2.5  —  4.0)  x  10*  to  accumulate  averages  to  calculate  the  density  profiles  n(  r) 
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and  the  mean  electrostatic  potential  0(r).  The  density  n(z)  was  obtained  by  counting  the 
ntimber  of  particles  in  slices  parallel  to  the  walls  and  dividing  by  the  volume.  Very  long 
nms  are  reqiiired  compared  to  the  primitive  electrolytes  simulations,  in  order  to  eliminate 
the  noise  in  the  density  profiles,  ^(r)  is  obtained  vta  the  relation 


“  T  /  dzi{z~zx  )qjnj(^zi ). 


Here  we  take  the  right  wall  as  the  reference  point  of  the  potential. 


(8) 


rV.  Results  and  Discussion 


In  the  following  we  report  our  simulation  results  of  the  density  profiles  and  mean 
electrostatic  potentials  for  the  different  systems  that  are  listed  in  table  I.  For  all  the 
systems  studied,  the  left  side  wall  is  charged  and  the  right  wall  (at  z  =$  )  is  neutral.  The 

surface  charge  density,  the  ion  and  solvent  density,  and  the  mean  electrostatic  potential 
are  reported  in  the  dimensionless  forms 

,  n,-  =  — ,  and  0*  =  (9) 

e  Hio 

where  n,o  is  the  bulk  number  density  for  species  i  ion  and  is  the  niunber  density  for  tne 
neutral  solvent. 

To  compare  the  results  obtained  here  with  those  obtained  for  the  primitive  electrdytes, 
in  figure  2,  we  show  the  SPM  resxilts  as  continuous,  dotted  or  dashed  curves  and  the  PM 
restilts  as  open  and  filled  circles  for  the  case  of  a*  si  0.42  and  C  =  IM.  These  PM  results 
were  obtsdned  in  step  (2)  of  the  simulations  and  agree  with  previotis  simulations  [4,7].  As 
is  well  known,  the  Gouy-Chapman  theory  and  simulations  of  the  primitive  electrolytes  are 
in  agreement  for  low  surface  charge  density  and  concentration.  Both  the  density  profile 
and  the  mean  electrostatic  potential  are  monotonic.  FVom  the  figtire  we  can  see  that,  the 
solvent  structure  included  in  the  SPM  simulation  has  a  dramatic  effect  on  the  electrolyte 
ion  profiles.  The  presence  of  solvent  molecules  at  the  walls  induces  strong  structure  in 
the  ion  distributions  at  both  charged  and  neutral  walls.  Near  the  left  wall,  the  neutral 
particles  and  counterions  have  five  layers,  the  coions  show  four  layers  and  near  the  right 
wall  there  are  four  layers  for  all  particles.  A  small  residue  of  the  coions  near  the  left  wall  is 
in  the  SPM  model  simulation.  These  structural  features  are  absent  in  the  PM  simtilations. 
Figure  3  shows  the  mean  electrostatic  potential  for  the  same  system  as  in  figure  2,  the 
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potential  decreases  rapidly  near  the  charged  wall  and  decays  to  zero.  As  mentioned  in 
section  III,  the  noise  in  the  density  profiles  prevents  us  from  getting  more  accurate  results 
to  see  a  clear  trend  of  the  potential. 

Figures  4  and  5  are  from  the  simulations  with  o*  =  0.7  and  C  =  IM.  The  layering 
of  ions  is  very  distinct  as  seen  in  figure  2.  Because  the  strong  field  on  the  left  wall,  the 
colons  were  completely  repelled  from  the  wall  in  the  PM.  However,  the  hard  sphere  solvent 
and  ion  interactions  in  the  SPM  again  result  in  a  significant  ordering  of  colons  near  the 
charged  surface.  The  mean  electrostatic  potential  of  the  SPM  decreases  more  rapidly  near 
the  charged  wall  than  does  that  of  the  PM.  The  deviations  of  our  simulations  are  shown 
in  errorbars  in  both  figures  4  and  5.  This  was  obtained  from  five  equal  length  segments 
of  the  simulations  with  different  seeds  for  the  random  number  generate  and  started  with 
the  previous  final  configuration  as  the  initial  configuration  of  the  new  run.  Notice  that  the 
errorbars  for  the  counterions  is  O.ld  ahead  in  its  z  coordinates  than  the  colons  and  solvent 
particles  in  figure  4. 

In  figures  6  and  7,  we  change  the  ion  concentration  to  2M  and  use  tr*  ^  0.4.  Again  the 
layering  of  ions  is  very  strong  and  similar  to  figures  2  and  3.  Furthermore,  a  charge  inversion 
occured  in  the  second  and  third  layer  with  the  colon  density  exceeding  the  covinterion. 

To  see  intuitively  how  a  dense  solvent  will  induce  significant  structure  in  the  ion 
density  distributions  even  thoxigh  the  ions  are  present  at  low  concentration,  let  us  consider 
a  dense  hard  sphere  mixture  confined  by  a  hard  wall.  Suppose  one  component  is  present 
at  high  density  and  the  others  are  present  at  low  concentrations.  Suppose  also  that  all 
hard  sphere  species  have  the  same  diameter.  Then  the  total  density  distribution  n(z) 
(n  =  ria)  near  the  wall  would  be  that  of  a  one  component  hard  sphere  fitiid.  At  high 
densities  this  is  known  to  be  a  highly  oscillatory  function.  The  probability  density  that 
any  hard  sphere  will  be  found  at  position  z  is  n{z)/N  (N  =s  Na)  and,  since  the  hard 
spheres  all  have  the  same  diameter,  the  probability  that  this  hard  sphere  is  species  a  is 
Xa,  where  Xa  is  the  mole  fraction  of  a.  Thus,  the  probable  density  na(r)  of  species  a  at 
z  is  na(z)  =  Nxain(z)/N)  =  Xon(z).  This  result  shows  that  the  density  profiles  of  all 
hard  sphere  species  have  the  same  structure,  the  difference  being  simply  that  their  ratios 
na(z)/nji(z)  scale  as  their  mole  fraction  ratio  Xa/x^j.  As  indicated  by  the  results  shown 
in  Figures  3-7,  adding  charge  to  the  confined  waUs  and  to  the  electrol3rte  modulates  the 
density  oscillations  somewhat,  but  does  not  appreciably  shift  the  positions  of  the  maxima 
and  minima  from  those  induced  by  the  hard  sphere  repulsions. 
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To  examine  the  effects  of  ion  diameters  on  the  highly  structured  ion  distributions  of 
the  SPM,  we  ran  simvilations  for  two  different  electrolytes  comprised  of  dissimilar  sized 
particles,  namely,  =  0.75d_  and  d+  =  0.5d_  .  The  results  are  given  in  figures  8 
through  11.  These  restilts  show  that  the  coion  peaks  shift  due  to  the  size  effect  but  there 
is  still  significant  layering  structure  in  the  ion  and  solvent  density  distributions.  Ion-size 
induced  charge  separations  is  cleiurly  seen  near  the  neutral  wall.  Due  to  the  noise  in  the 
density  profiles,  we  cannot  msdce  an  accurate  estimate  of  the  potentials  resulting  from 
these  separations.  Comparing  the  density  profiles  in  figures  8  and  10,  one  can  see  that 
the  layering  seems  mainly  determined  by  the  larger  solvent  and  covmterions  near  the  two 
walls,  even  though  the  effect  of  the  coion  size  made  the  peak  shift  a  little.  Again  these 
structures  are  absent  in  PM  electrolytes. 

From  the  above  analysis  and  comparision  with  the  primitive  electrolyte  simulations 
we  see  the  essential  role  of  the  solvent  in  determining  the  structme  of  electrolytes  at  solid 
surfaces.  In  fact,  the  layering  structures  appearing  in  the  SPM  electrolytes  are  mainly  due 
to  the  high  total  fluid  density  resulting  from  the  presence  of  the  solvent  molecules.  The 
small  residue  of  coions  near  the  charged  surfaces  is  due  to  the  layering  of  the  solvent  and 
counterions.  The  fact  that  we  did  not  observe  stronger  structures  near  the  charged  wall  at 
higher  surface  charge  density  (<r*  =  0.7)  than  at  lower  surface  charge  densities  ( <7*  ~  0.4) 
confirms  the  idea  that  the  effect  of  the  solvent  plays  a  more  important  role  in  determining 
the  electrolyte  structure  than  solely  by  the  surface  charge  hard  sphere  ion  interactions,  as 
is  the  case  of  the  primitive  electrolyte.  In  all  cases  that  have  been  studied  in  this  paper, 
the  mean  electrostatic  potential  of  the  SPM  is  smaller  than  that  obtained  &om  the  PM 
near  the  charged  wall. 

Finally,  we  test  our  results  for  consistency  by  comparing  in  Table  II  the  simulated 
difference,  0(0)  —  0(d/2),  with  the  value  of  2s’<7£i/«,  since  theoretically 

0(0)  =  0(d/2)  +  ^  .  (10) 

The  numbers  in  Table  II  show  that  the  agreement  is  very  good.  Also  shown  in  Table  II 
are  values  obtained  by  Valleau  and  coworkers  in  a  grand  canonical  ensemble  simulation  of 
the  primitive  model  electrolyte  at  the  same  conditions  that  we  used.  Their  restilts  also 
obey  eq.(lO)  quite  accurately.  Our  values  of  0(0)  and  0(d/2)  differ  somewhat  from  theirs. 
The  difference  perhaps  arises  from  size  difference  of  the  simulated  systems  (they  use  much 
larger  systems  than  ours). 
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The  input  surface  charge  density  tr  agreed  with  that  computed  from  the  ion  density 
distributions, 

I  tfi 

<7  =  —/  y'9<n,(r)dz,  (11) 

A  Jdit  ^ 

to  the  fourth  digit.  This  is  guaranteed  in  our  simulations  by  the  charge  balance  of  eq.(l}. 

The  apparent  structure  of  the  interfaces  obtained  from  the  solvent  primitive  model 
may  provide  new  insight  into  the  interaction  of  charged  particles  near  siufaces.  The  layering 
of  molecules  near  a  surface  would  be  reflected  by  oscillatory  interaction  forces.  Force 
measurements  using  mica  surfaces  have  revealed  short  range  oscillatory  electrical  double 
layer  forces  [25-29]. 
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Table  I.  Parameters  used  in  this  paper,  =  a^fc,  side  length  £,  surface  separation 
H.  Number  of  cations,  anions,  and  neutral  particles  are  ,  iV_,  and  No,  and  Nc  — 
N+  +  iV_,  AN  =  N-  —  iV+  =  erA,  A  =  L  x  L.  C  is  the  bulk  ion  concentration  for  PM 
electrolytes  in  molar.  The  number  of  steps  run  in  the  simulations  are  given  as  Nsteps  in 
unit  of  10®  steps  in  the  last  column.  In  all  the  cases  studied  in  this  paper  the  neutral 
solvent  particle  and  the  zmion  have  diameter  of  d©  =  d-  =  d  =  4.25  A.  For  systems  a-6,  c, 
d  and  c,  ln7±  =  —0.127,  0.271,  —0.318  and  —0.44,  respectively. 


sys 

H/d 

L/d 

d-^/d- 

Nc 

No 

AN 

C(M) 

Nsteps/ 10^ 

a 

12.0 

0.42 

26 

134 

8 

1.0 

4.47 

0.70 

32 

136 

14 

1.0 

2.5 

4.5 

1.0 

0.40 

54 

146 

8 

2.0 

2.6 

12.0 

4.36 

0.75 

0.42 

26 

134 

8 

1.0 

2.5 

12.0 

4.36 

0.5 

0.42 

26 

134 

8 

1.0 

2.5 

Table  II.  Results  of  the  SPM  (computed  in  canonical  ensemble,  CMC)  and  PM  (com¬ 
puted  in  the  grand  canonical  ensemble, GCMC)  simulations  of  the  total  potential  drop 
across  the  double  layer  V’o  =  ^**(0)  diffuse  layer  mean  electrostatic  potential 

i/?*(d/2)  as  calculated  from  eq(8)  with  a  bin  size  dz  =  O.ld  and  compared  with  eq(ll).  For 
systems  d  and  e,  we  give  the  diffuse  layer  potential  at  d/2.  The  numbers  in  parenthesis 
sure  the  numerical' uncertainties  determined  from  consequent  simulation  segments. 


sys 

a 

b 

c 

d 

e 

SPM  0*(f) 

1.11 

1.75(0.18) 

0.45 

1.06 

1.32 

00 

5.52 

9.10(0.18) 

4.60 

5.47 

5.73 

r^-rd) 

4.41 

7.35 

4.15 

4.41 

4.41 

PM  0*(|) 

2.78(0.04) 

4.60(0.02) 

1.75(0.07) 

2.82(0.02) 

2.89(0.02) 

3.08(0. l)t 

5.7l(0.14)t 

2.29(0.09)t 

05 

7.19(0.13) 

11.95(0.04) 

5.90(0.07) 

7.23(0.02) 

7.30(0.02) 

7.52t 

13.10t 

6.47t 

«  -  r(j) 

4.41 

7.35 

4.15 

4.41 

4.41 

4.44t 

7.39t 

4.18t 

litari 

_ t _ 

4.41 

7.356 

4.16 

4.41 

4.41 

t  GCMC  simtilations  from  reference  [5], 
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Figure  Caption* 


Figure  1. 
Figtu-e  2. 

Figure  3. 
Figure  4. 

Figure  5. 
Figxire  6. 
Figure  7. 
Figure  8. 
Figure  9. 
Figure  10. 
Figure  11. 


Schematic  representation  of  solvent  moiectdes  and  ions  of  diameter  d  near  a  uniformly 
charged  flat,  impenetrable  wall. 

Reduced  Density  profiles  n*(z)  =  n{z)/no  for  a  1:1  electrolyte  at  C=lM  and  d+  = 
d-  =  d  s=  4.25  A.  The  left  wall  has  a  sruface  charge  density  of  a*  =  ft  —  0.42  and 
the  ri^t  wall  is  unchrirged. 

Reduce  Mean  electrostatic  potential  ^*{z)  for  a  1:1  electrolyte  for  the  same  conditions 
as  in  figure  1. 

Reduced  Density  profiles  n*(2)  for  a  1:1  electrolyte  at  C=1M  and  a*  =  0.7  .  Other 
conditions  emd  keys  are  the  same  as  in  figure  1.  The  errorbars  for  the  counterions  are 
O.ld  ahead  in  the  z  coordinates  than  the  colons  and  solvent  particles. 

Reduce  Mean  electrostatic  potential  ^*{z)  for  a  1:1  electrolyte  for  the  same  conditions 
£is  in  figure  3.  Errorbars  are  shown  for  SPM. 

Reduced  Density  profiles  n*{z)  for  a  1:1  electrolyte  at  Cs=2M  and  =  0.40  .  Other 
conditions  and  keys  are  the  same  as  in  figure  1. 

Reduce  Mean  electrostatic  potential  ^*{z)  for  a  1:1  electrolyte  for  the  same  conditions 
as  in  figxxre  5. 

Reduced  Density  profiles  n*{z)  for  a  1:1  electrolyte  at  Ca=lM  and  =  0.42  .  d+  = 
|d_.  Other  conditions  and  keys  are  the  same  as  in  figure  1. 

Reduce  Mean  electrostatic  potential  0*(r)  for  a  1:1  electrolyte  for  the  same  conditions 
as  in  figure  9. 

Reduced  Density  profiles  n*(z)  for  a  1:1  electrolyte  at  C=1M  and  a*  =  0.42  .  d+  = 
^d_.  Other  conditions  and  keys  are  the  same  as  in  figure  1. 

Reduce  Mean  electrostatic  potential  for  a  1:1  electrolyte  for  the  same  conditions 
as  in  figrire  7. 
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